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^ ' Abstract. As a part of an ongoing effort to characterise the young stellar populations in the Large Magellanic Cloud, we 

, present HST-WFPC2 broad and narrow band imaging of two fields with recent star formation activity in the Tarantula region. 

A population of objects with Ho- and/or Balmer continuum excess was identified. On account of the intense Ho- emission 
' - ^ ' (equivalent widths up to several tens of A), its correlation with the Balmer continuum excess and the stars' location on the HR 

. diagram, we interpret them as low mass (~ 1 - 2 Mq) Pre-Main Sequence stars. In this framework, the data show that coeval 

• high and low mass stars have significantly different spatial distributions, implying that star formation processes for different 

^— ■») I ranges of stellar masses are rather different and/or require different initial conditions. We find that the overall slope of the mass 

■ function of the young population is somewhat steeper than the classical Salpeter value and that the star formation density of this 
i-^ ' young component is 0.2 - 0.4 MQyr^'kpc"^, i.e. intermediate between the value for an active spiral disk and that of a starburst 

region. The uncertainties associated with the determination of the slope of the mass function and the star formation density are 
^ ■ thoroughly discussed. 

, Key words. Galaxies: individual (Large Magellanic Cloud) - Galaxies: evolution - Stars: fundamental parameters - Stars: 

' formation - Stars: pre-main sequence - Stars: mass function 

> : 
•I— I 

1 . Introduction Padget 1996). Also, they typically cover large areas on the sky 

jrt that require extensive tiling to be fully covered with the current 

. - - Star-forming processes determine much of the appearance of ^^^^^^^^^ ^^j^^ ^^^^^^ ^^^^ ^^^^-^^ ^^^-^^^ example, covers 
the visible Universe: the shape of the Initial Mass Function ^^^j^ ^2 degrees on the sky). 
(IMF) and its normalisation (the star-formation rate) are, to- 



gether with stellar evolution theory, key ingredients in deter 
mining the chemical evolution of a galaxy and its stellar con- 
tent. Yet, our theoretical understanding of the processes that 



However, studying the effects of lower metallicity on star 
formation is essential to understand the evolution of both our 
own Galaxy, in which a large fraction of stars were formed at 



1 J r ' j vcr 1 1 1 J . . • .11 . . metallicities below solar, and of what is observed at high red- 
lead from dirruse molecular clouds to stars is still very tenta- 

, , . , , • 1. 1 shifts. As a matter of fact, the global star formation rate ap- 

tive, as many complex physical phenomena are simultaneously uu u ■ if fir\ \ 



pears to have been much more vigorous (a factor of 10 or so) 

at z ^ 1.5 than it is today (Madau et al 1996 and subsequent 

incarnations of the so-called "Madau plot")- At this epoch the 

, , , . , , , , , , . mean metallicity of the interstellar gas was similar to that of the 

It, on the one hand, studying them leads to the the best possi- . „ . „, , ,r ..^ r. • r- n o n 1^,r,^,^ 

. . ... . Large Magellanic Cloud (LMC, e.g. Pet, Fall & Hauser 1999). 



at play in producing the final results. 

From an observational standpoint, most of the effort was 
traditionally devoted to nearby Galactic star-forming regions. 



ble angular resolution, on the other hand, this is achieved at the 
expense of probing only a very limited set of astrophysical con- 
ditions (all these clouds have essentially solar metallicity, e.g 



This fact combined with the galaxy's well-determined distance 
and small extinction makes the study of star forming regions in 
the LMC an important step towards understanding the general 
* Based on observations with the NASA/ESA Hubble Space picture of galaxy evolution. 
Telescope, obtained at the Space Telescope Science Institute, which With a distance modulus of 18.57 + 0.05 (see the discussion 
is operated by AURA, Inc., under NASA contract NAS 5-26555. in Romaniello et al 2000), the LMC is the closest galactic com- 
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panion to the Milky Way after the Sagittarius dwarf galaxy. At 
this distance one arcminute corresponds to roughly 15 pc and, 
thus, one pointing with most of the current generation of in- 
struments comfortably covers almost any individual star form- 
ing region in the LMC (see, e.g., Hodge 1988). In particular, 
the field of view of 2.5' x 2.5' of the camera we have used, 
the WFPC2 on board the HST, corresponds to 37 x 37 pc^ and 
leads to the detection of several thousands of stars per pointing. 
The LMC is especially suited for stellar populations studies for 
two additional reasons. First, its depth along the line of sight 
is negligible, at least in the central parts we consider (van der 
Marel & Cioni 2001), and all the stars can be effectively con- 
sidered at the same distance, thus minimising a spurious scatter 
in the Colour-Magnitude Diagrams (CMDs). Second, the ex- 
tinction in its direction due to dust in our Galaxy is low, about 
E(B - V) ^ 0.05 (Bessell 1991, Schwering 1991), and hence, 
our view is not severely obstructed. 

The advent of the Hubble Space Telescope (HST) has for 

the first time made it possible to observationally tackle the open 
questions about star formation in outer galaxies and, in partic- 
ular, in the LMC, down to a solar mass or even lower. Some 
of the evidence from ground and HST-based studies shows 
that there may be significant differences between star forma- 
tion processes in the LMC and in the Galaxy. Immediately af- 
ter the first HST refurbishment mission, the observations of 
the double cluster NGC 1850 in the LMC made with Wide 
Field Camera 2 (WFPC2) have provided the first detection of 
a population of Pre-Main-Sequence (PMS) stars (Gilmozzi et 
al 1994) in an extragalactic star forming region. The evidence 
was purely statistical, in the sense that the existence of PMS 
stars was deduced by the presence of many stars that were 
lying above the main sequence and whose number (almost 
400) could not be accounted for by any evolved population. 
Subsequently, Romaniello (1998), Panagia et al. (2000) and 
Romaniello et al (2002), by comparison of the magnitudes ob- 
tained with an Ha narrow band filter (F656N) with those taken 
with a broad band red filter(F675W), measured the equivalent 
width of the Ha emission for several hundred of candidates 
PMS stars in the field around SN 1987 A, thus confirming their 
PMS identification unambiguously. Eventually, the first spec- 
troscopically confirmed discovery of a bona-fide T Tauri star 
in the LMC, LTS J054427-692659, a low-mass, late-type star 
located within the dark cloud Hodge II 139, was reported by 
Wichmann, Schmitt & Krautter (2001). 

Some of the evidence from ground and HST-based studies 
shows that there may be significant differences between star 
formation processes in the LMC and in the Galaxy. For exam- 
ple, Lamers et al (1999) and de Wit et al (2002), on the basis of 
ground-based data, have suggested the presence of high-mass 
Pre-Main Sequence stars in the LMC (Herbig AeBe stars, but 
see the caveats in de Wit et al 2005) with luminosities system- 
atically higher than observed in our Galaxy, and located well 
above the "birthline" of Palla & Stabler (I99I). They attribute 
this finding either to a shorter accretion timescale in the LMC 
or to its smaller dust-to-gas ratio. Whether such differences in 
the physical conditions under which stars form will generally 
lead to differences at the low mass end is an open question, but 



Panagia et al (2000, Paper I) and Romaniello et al (2004) offer 
tantaUsing evidence in this direction. 

HST-based studies of the IMF in the LMC (which have 
mostly concentrated on young, compact clusters) have in fact 
produced, for the lower masses, often discrepant results. On the 
one hand, there is widespread agreement that the IMF for M > 
3Mq is similar to the mean Galactic one (i.e. a "Salpeter func- 
tion", a power-law with slope a = -2.35). On the other hand, 
though, the results for the lower masses (IMq < M < 3Mq) 
are rather difficult to reconcile with each other, with different 
authors finding (in different regions, or even in the same one) 
very different slopes, ranging from very steep IMFs (a ^ -3, 
Mateo 1988) to quite shallow ones (a ^ - 1 , Sirianni et al 2000) 
with many intermediate values represented. 

As a part of an ongoing effort to characterise the young 
stellar population in the LMC, we have applied the techniques 
developed by Romaniello (1998) to detect and characterise low 
mass Pre-Main Sequence (PMS) stars in two regions of recent 
star formation observed with the WFPC2 on board the HST. 
The fields were selected in the HST archive so as to have deep 
narrow-band Ha observations, which, as we shall see, are fun- 
damental to detect the low mass Pre-Main Sequence popula- 
tions. Both fields are in the proximity of Supernova 1987A 
(SN1987A), whose immediate vicinities we have already stud- 
ied with similar methods (Romaniello 1998, Panagia et al 2000, 
Romaniello et al 2004). In addition to Ha, one of the two fields 
presented here was imaged in four broad band filters, so that 
this wealth of data allowed us to recover the intrinsic proper- 
ties of the stars as well as their reddening following the pre- 
scriptions developed by Romaniello et al (2002). For the other 
field the broad band data were limited to two filters only, and, 
therefore, the analysis had to be necessarily cruder. We will, 
however, show how interesting conclusions on the properties 
of the stellar populations can be drawn also from such a lim- 
ited dataset. More fields are in the process of being analysed 
to enlarge our study of young populations in the Magellanic 
Clouds across the mass spectrum. 

Section 2 is devoted to the description of the datasets and 
the reduction techniques, whereas the results for the two fields 
considered here are discussed in section 3 and 4. Finally, sum- 
mary and conclusions are presented in section 5. 

2. Observations and data reduction 

We have taken advantage of the vast HST archive to select two 
LMC fields with recent star formation activity. They were cho- 
sen for having both broad-band and Ha observations, which, 
as we shall see in section 3.4, is a fundamental requirement to 
identify the low mass Pre-Main Sequence stars, and unambigu- 
ously disentangle them from the much older field subgiants. 

The two fields we present in this paper are located some 20' 
southwest of the 30 Doradus nebula. This area includes regions 
of very active star formation, in which different groups of early 
type stars are interspersed with HII regions and Supernova rem- 
nant shells, as well as one of the best know objects in the entire 
sky: Supernova 1987A. The footprints of the WFPC2 point- 
ings superimposed on a DSS image are shown in Figure 1. 
Also shown in the same Figure are the locations of the fields 
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around SN1987A itself in which Panagia et al (2000) and and 
Romaniello et al (2004) have identified nearly 500 Pre-Main 
Sequence candidates through their Ha and Balmer continuum 
emission. 



Filter Name 


Exposure Time (s) 


Comments 


F547M 


2 X 120" 


V-like 


F675W 


2 X 120" 


R-like 


F656N 


1000 + 2 X 2000" 


Ha 




900 + 1000 + 2000" 






Fig. 1. Footprints of the WFPC2 pointings superposed to a DSS image 
20' on the side. The progenitor of SN1987A is visible at the centre 
of the three WFPC2 fields of Panagia et al (2000). The 30 Doradus 
nebula is roughly 20' to the northeast. 



The centres of the fields are located at a = 05 : 36 : 1 1.57, 
(5 = -69 : 23 : 00.74 (J2000) and a- = 05 : 34 : 38.01, 
6 - -69 : 17 : 14.90 (J2000). Given their position with respect 
to the SN1987A, in the following we will refer to the pointings 
described here as South field and West field, respectively. Their 
angular separations from SN1987A are 8' and 4.5' (i.e. 120 
and 68 pc in projection), respectively, and the South field is 
just 2.5' west of the young cluster NGC 2050. The logs of the 
observations of the two fields are reported in Table 1 and 2. The 
detailed description of the WFPC2 camera and its filters can 
be found in the corresponding Instrument Handbook (Heyer et 
al 2004). 



Filter Name 


Exposure Time (s) 


Comments 


F300W 


600 + 1200 


WideU 


F4S0W 


50 + 200 


WideB 


F675W 


50 + 200 


R-like 


F814W 


140 + 400 


l-like 


F656N 


1200 + 1200 + 1300 


Ha 



Table 1. Log of the observations of the South field. They were taken 
on January, P' 1997, under proposal number 6437, PI Robert P. 
Kirshner. 



The data were processed through the standard Post 
Observation Data Processing System pipeline for bias removal 



" June, 24"> 1996. 
* June, 25*1996. 

Table 2. Log of the observations of the West field. The proposal num- 
ber is 6033, PI Jeremy Walsh. 



and flat fielding. In all cases the cosmic ray events were re- 
moved combining the available images, duly registered. 

The plate scale of the camera is 0.045 and 0.099 arc- 
sec/pixel in the PC and in the three WF chips, respectively. 
We performed aperture photometry following the prescriptions 
by Gilmozzi (1990) as refined by Romaniello (1998), i.e. mea- 
suring the flux in a circular aperture of 2 pixels radius and 
the sky background value in an annulus of internal radius 
3 pixels and width 2 pixels. Due to the undersampling of the 
WFPC2 PSF, this prescription leads to a smaller dispersion in 
the Colour-Magnitude Diagrams, i.e. better photometry, than 
PSF fitting for non-jittered observations of marginally crowded 
fields (Cool & King 1995, Romaniello 1998). The photometric 
error is computed taking into account the poissonian statistics 
of the flux from the source and the measured rms of the flux 
from the sky background around it. This latter includes the ef- 
fects of poissonian fluctuations on the background emission, 
flat-fielding errors, etc. 

Photometry for the saturated stars was recovered by either 
fitting the unsaturated wings of the PSF for stars with no sat- 
uration outside the central 2 pixel radius, or by following the 
method developed by GilUland (1994) for the heavily saturated 
ones. 

The flux calibration was obtained using the internal caUbra- 
tion of the WFPC2 (Whitmore 1995), which is typically accu- 
rate to within 5% at optical wavelengths. The spectrum of Vega 
is used to set the photometric zeropoints (VEGAMAG system). 

3. The South field 

The F450W (B), F675W (R) and F656N (Hor) composite image 
of the South field is shown in Figure 2. 

We have detected 13,098 stars in the whole field, 4,108 
of which with a mean error in the 4 broad bands lower than 
0.1 mag: 



54 



6^ 



'V45()W ^ 'VhVSW 



- (5r- 



< 0.1 



(1) 



The brightest star in the field, located in the WF3 chip, is 
so badly saturated that it was impossible to measure its mag- 
nitude from the WFPC2 images. On the other hand, it is a 
well known star, Sk -69 211 (Sanduleak 1969), also known 
as HDE 269832, and its accurate photometry is available from 
ground-based studies. Fitzpatrick (1988) classifies it as B9Ia 
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Fig. 2. The South field as observed in the combination of F450W (B), 
F675W (R) and F656N (Ha) filters. 



supergiant with V= 10.36, B-V= 0.09 and E(B - V) ^ 0.1. 
With a radial velocity of +272 km s~' it is a bona fide mem- 
ber of the LMC (Ardeberg et al 1972). Using the temperature 
scale and bolometric corrections of Schmidt-Kaler (1982), one 
derives for it an effective temperature of roughly 10,000 K and 
a luminosity of roughly 3-10^ L0. 

The Colour-Magnitude Diagrams for four combination of 
filters are shown in Figure 3. 

The presence of at least two distinct populations is read- 
ily detected in Figure 3. The first one is associated with 
the feature extending up to mF45ow - 15 and with a colour 
(mF45ow - 'wpSMw) - -0.6 (see panel (c)). Of course, this is a 
young generation of stars still on or just off the Zero Age Main- 
Sequence. The old population of the LMC field, with ages in 
excess of 500 Myr or so, is identified by the presence of the 
Red Clump at otf45ow - 20, (mF45ow - ffiF8i4w) - 1-8 and an 
extended Red Giant Branch reaching magnitudes as bright as 

'WF450W - 18 at (OTF450W - 'WF814w) - 3. 

3.1. Reddening distribution 

In order to derive the detailed properties of the stars in the 
field we first have to determine the interstellar reddening. As 
in any region of recent star formation, one can expect the stars 
to be still entwined with the material they were formed from 
and, hence, the interstellar reddening to be spatially inhomo- 
geneous. We have, then, applied the procedure developed by 
Romaniello et al (2002) to measure E(B - V), T^// and lumi- 
nosities for individual stars. We have used a value for the dis- 
tance to the LMC of 5 1 .8 kpc (Romaniello et al 2000). In brief, 
the procedure consists in a fit of the observed fluxes to a grid of 
theoretical ones, reddened by various amounts of E(B - V). A 




Fig. 3. Colour-Magnitude diagrams of the 13,098 stars detected in the 
South field for four combination of filters. The black dots are the 4,108 

stars with a mean photometric error in the four bands (^4, see equa- 
tion 1) smaller than 0.1 mag. Sk -69 211 is not shown. 



minimum technique is, then, used to identify the best combi- 
nation of Teff, E(B - V) and luminosities, with their associated 
errors. As inputs for the fit we have used the atmosphere mod- 
els by Bessell et al (1998) and the reddening law appropriate 
for this region of the LMC as derived by Scuderi et al (1996) 
from a study of Star 2, one of the two stars projected within a 
few arcseconds from SN1987A. 

We were able to estimate E(B - V) for 1,205 stars, i.e. 9% 
of the total or one every 4 square arcseconds on average, while 
for the others we have used the mean value of their neigh- 
bours with direct measurements (cf. Romaniello et al 2002). 
An inspection of Figure 10, which will be discussed in detail 
in section 3.3.2, shows that the location of the dereddened stars 
in the (mpsoow.o - 'Wf450w,o) vs (otf450w,o - wpgww.o) colour- 
colour plane generally agrees very well with the expectations 
from the BesseU et al (1998) model atmospheres, with a scatter 
consistent with observational errors (^4 < 0.1, the discrepancy 
observed for some stars at (otf450w,o - »if814w,o) - 0.7 is thor- 
oughly discussed in section 3.3.2). This agreement, then, gives 
us confidence on the use of the Bessel et al (1998) atmospheres 
to reproduce the observed colors and confirms that our dered- 
dening procedure does not introduce any significant systematic 
effect (see also Romaniello et al 2002). 

The derived reddening distribution is indeed rather clumpy, 
as shown in Figure 4. The mean extinction over the field is 
E(B - V) = 0.180 (second column of Table 3), i.e. compa- 
rable to the mode of the distribution for the old population 
of the LMC as a whole of E(B - V) ^ 0.1 (Zaritsky 1999). 
The corresponding rms is 0.091 (third column). Thanks to the 
large number of stars, then, the error on the mean is negli- 
gible, i.e. 0.091/ V(1205) = 0.003. Again from Table 3, one 
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sees that the mean formal error on the individual reddening 
determination due to the combined effects of photometric er- 
rors and fitting uncertainties (fourth column) is 0.029 magni- 
tudes. This, together with the measured rms, implies that the 
intrinsic dispersion in E(B - V) due to the patchy distribu- 
tion of dust is 0.086 magnitudes. This value is both statisti- 
cally significant and non negligible, considering that it reflects 
itself in a spread of about 0.086 x 5.5 0.5 mag in F300W, 
0.086 X 3.9 - 0.3 mag in F450W, 0.086 x 2.4 ^ 0.2 mag 
F675W and 0.086 x 1.9 ^ 0.15 mag in F814W (see Table 8 
of Romaniello et al 2002 for the extinction coefficients in the 
WFPC2 bands). 



Chip Mean rms Error Dispersion No. of stars 

0233 OTn 0037 OTTI 33 

WF2 0.210 O093 O027 0.089 401 

WF3 0.157 O089 O031 0.083 427 

WF4 0.168 0.077 0.028 0.072 344 

All O180 0.091 0.029 0.086 U05 

Table 3. Spatial distribution of E(B - V) colour excess as measured 
in the South field. In column 3 we report the measured rms of the 
E(B - V) distribution, in column 4 the mean error on E(B - V) result- 
ing from the fitting procedure (see text), whereas column 5 lists the 
intrinsic reddening dispersion, computed as the quadratic difference 
of the first two quantities. The number of stars with direct extinction 
determination is shown in the last column. 
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Fig. 4. Spatial distribution of stars with high (E(B - V) > 0.3, filled 
circles) and low (E(B - V) < 0.1, open circles) reddening. The av- 
erage reddening is significantly higher in the WF2 chip, i.e. towards 
NGC 2050 (see Figure 1). 



3.2. HR diagram 

Once dereddened, the stars can accurately be placed in the HR 
diagram and their properties derived by comparing their posi- 
tions with theoretical models (Figure 5). 
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It is also clear from Figure 4 that the interstellar extinction 
is inhomogeneous on larger scales as well, with the reddening 
being significantly higher in the PC and WF2 chips, i.e. in the 
direction of NGC 2050, than elsewhere. The actual numbers 
are reported in Table 3, where one can see that the difl'erence in 
the mean reddening between the PC and WF2 chips on the one 
side and the WF3 and WF4 chips on the other is statistically 
highly significant. The dispersion, on the other hand, is almost 
constant across the chips (see also section 3.5.1). All of these 
facts highUght once more how crucial it is to accurately dered- 
den the stars in order to recover and interpret their properties. 

Finally, it is interesting to notice that the spatial distribution 
of E(B - V) does not correlate with the difliise Ha emission, 
which is highest across the PC and the upper part of the WF3 
chips, whereas it is quite faint in the WF2 chip (see Figure 2). 



Fig. 5. Hess diagram for the 13,098 stars in the South field. The open 
squares indicate 6 stars for which the F300W magnitude is ill deter- 
mined because of saturation and, hence, the fit was performed exclud- 
ing this filter. All of them, but one, require a temperature equal to 
or larger than 50,000 K, the highest available for the model atmo- 
spheres by Bessell et al (1998). The location of Sk -69 211 according 
to the photometry by Fitzpatrick (1988) is shown with a star sym- 
bol. The ZAMS and selected isochrones for Z= 0.3 ■ Zq by Brocato & 
Castellani (1993) and Cassisi et al (1993) are overlaid to the data (grey 
soUd lines). These models only extend to 25 Mq, so we have comple- 
mented them with those of Schaerer et al (1993) for higher masses 
(grey dashed lines). 

As we have already noticed from an inspection of the 
CMDs shown in Figure 3, no single age can explain the dis- 
tribution of stars as observed in the HR diagram of the South 



6 



Martino Romaniello et al.: Pre-Main Sequence stars in the LMC 



field. Starting from the top-left corner of the HR diagram, the 
most massive stars, plotted as open squares in Figure 5, have 
masses of the order of 40 Mq and, thus, are definitely younger 
than 5 Myr. Since they are severely affected by saturation in 
both F300W exposures, their intrinsic parameters were derived 
excluding the flux in this filter from the fit. As discussed by 
Romaniello et al (2002), lacking the UV information, the fit is 
rather uncertain. In any case, all of them, but one, seem to re- 
quire temperatures higher than 50,000 K, the highest available 
for the Bessell et al (1998) models. In conclusion, although 
their stellar parameters cannot be accurately measured, these 
stars are certainly massive and, hence, young, even more so 
than their position in the HR diagram would indicate. 

With a comparable bolometric luminosity, but a much 
lower temperature, there is Sk -69 211. It falls exactly on the 
7 Myr isochrone by Schaerer et al (1993), i.e. it is without a 
doubt older than the bluest stars in the field. A comparison with 
stellar evolutionary models assigns a mass of roughly 27 Mq to 
it. Interestingly, a similar phenomenon was reported in Galactic 
OB associations by Massey et al (1995), who noted the "occa- 
sional presence of an evolved star" among a younger popula- 
tion. 

Finally, stars with ages between a few tens of million and 
several billion years are required to account for the stars at 
lower luminosities. 

3.3. A population of peculiar objects 

A population of peculiar objects was detected in the South field 
by means of their Ha and Balmer continuum excesses (the two 
emissions positively correlate, as we shall see in section 3.3.2). 
In this section we describe the selection criteria for these stars 
and offer a possible explanation of their nature as Pre-Main 
Sequence stars. 

3.3.1. Ha emission 

Following Romaniello (1998; see also Panagia et al 2000 and 
Romaniello et al 2004), we compare the magnitudes in the R 
band (F675W) with the ones measured with the narrow band 
Ho; filter (F656N) to identify stars with a sufficiently strong 
Ha emission. As it can be seen in Figure 2 there is significant 
nebular Ha emission in the South field. It is distributed in fila- 
ments with a typical length of several tens of arcseconds up to 
several arcminutes. All of the stars that are found to have Ha 
emission (see below) were inspected by eye to ensure that the 
nebular contribution was properly subtracted from the narrow 
band photometry. 

In order to select stars with a statistically significant 
Ha emission we have compared the observed histogram of 
(mp675w - »Jf656n) for the stars in the field with the distribution 
expected if no emission was present. This latter was derived 
with the following four steps: 

1. Kurucz (1993) model atmospheres for a metallicity 
[M/H] = -0.5 at different effective temperatures were used 
to compute synthetic spectra between 5,700 and 7,800 A, 
the region covered by the F675W filter. A spectral reso- 
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Fig. 6. Relation between the WFPC2 (mFeTSw-wpeseN) colour {i.e. R- 
Hor) and effective temperature as derived from Kurucz (1993) model 
atmospheres (see text). 



lution of 1 A was chosen, so as to sample the Ha line to 
a much better accuracy than the FWHM of the WFPC2 
F656N filter (~ 20 A, Heyer et al 2004). 

2. These spectra were fed to the synphot task in IRAF to 
yield the expected (mpevsw - 'wpfeseN) colour of "normal" 
stars, i.e. with Ha in absorption, as a function of their 
effective temperature. The resulting relation is shown in 
Figure 6. 

3. We have used this relation together with the measure effec- 
tive temperatures of the stars (cfr section 3.2) to compute 
the expected (mpevsw - wif656n) colour for our stars in the 
absence of any Ha emission (shaded area in Figure 7). 

4. Finally, this intrinsic distribution was broadened by the 
measured photometric error on the (;7tF675w - 'wf656n) 
colour. To this end, the intrinsic distribution was convolved 
with the appropriate error function, i.e. the normalised sum 
of the error Gaussians of the individual stars. This is, then, 
the colour distribution that the stars in the field are expected 
to have according to their measured temperature. 

The comparison of the observed and expected distributions 
derived as described above is presented in Figure 7 (soUd and 
dashed lines, respectively). The left-hand tails of the two dis- 
tributions are in good agreement with each other, implying that 
the broadening towards negative values of (?mf675w - 'WFesew) 
compared to the intrinsic, error-free distribution (shaded area 
in Figure 7) is caused by photometric errors. Moving towards 
more positive (otf675w - ^MpeseN) values, however, the two his- 
tograms become increasingly discrepant, until an excess of ob- 
served stars with positive (mpeTsw - wpeseN) colours, i.e. Ha 
emission, compared to the expectations is evident for (wf675w - 
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wf656n) ^ 0.2 (dots in Figure 7). Applying ax^ test confirms 
that the distribution are conclusively different (x^ = 11). 



600 



400 



3 
2 



200 




-1 -0.5 0.5 1 

m(F675W) - m(F656N) 

Fig. 7. Observed (thick solid line) and expected (thin dashed line, see 
text) distribution of (mpeTsw - "tpeseN) colour in the South field. The 
histograms contain 4,938 stars, those with a positive detection in the 
F656N (Ha). The shaded areas marks the intrinsic colour interval pre- 
dicted by the atmosphere models, i.e. without the effect of observa- 
tional errors (see also Figure 6). As it can be seen, the left-hand tails 
of the two histograms are in agreement, and the excess of observed 
stars with positive (/Mf675w - WpeseN) colours, i.e. Ha emission, com- 
pared to the expectations is evident for (mpeTsw - 'Mf656n) ^ 0.2. The 
distribution of the stars with significant excess (see text) is shown as 
dots, together with the corresponding errorbars. The value of the ex- 
cess above which stars can be identified individually is marked with a 
vertical line. 

In Figure 8 we report the observed and expected (wf675w - 
"JpeseN) distributions, derived as described above, broken down 
in different intervals of observational uncertainty dimpeisw - 
'Wf656n)- The error range used is quoted in each panel. Since 
the photometric error is dominated by Poisson noise on pho- 
ton statistics, the intervals in error are essentially intervals in 
^Mpevsw magnitude (reported within brackets in each panel). In 
all cases the widths of the left-hand side of the observed and 
expected distributions agree, confirming that the broadening is, 
indeed, due to photometric errors, whereas an excess of Ha- 
emitting stars is detected at (mpevsw - wf656n) ^ 0.2. 

Based on the distributions of Figure 7 and 8, we use the 
following criteria to identify the stars with a statistically signif- 
icant excess: 

- the stars that populate (mpeTsw - fnFesen) bins in which 
the observed distribution exceeds the expected one by 3 
Poissonian standard deviations or more; 

- while the criterion above ensures a high statistical signifi- 
cance of the selected objects, by definition it excludes stars 
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Fig. 8. Same as Figure 7, but each panel only contains the stars in the 

quoted range of error on the mpevsw - 'Mf656n colour. The correspond- 
ing approximate range in mpeTsw magnitude is also reported in each 
panel. The value of the excess above which stars can be identified in- 
dividually is marked with a vertical line in each panel. 

in scarcely populated bins, i.e. with less than 9 stars. To 
compensate for this, we also include in the sample the stars 
for which the (rnpevsw - otf656n) excess itself is highly sig- 
nificant, i.e. (WF675W - Wf656n) > 3 • (5(mF675W - Wf656n), 

even if they belong to sparse bins. This only adds roughly 
1% of the stars to the sample, but they are important ones 
because they have the highest Ha emission and, hence, are 
the most active ones. 

In total, there are 781 + 33 stars that satisfy the crite- 
ria above. Their distribution, with the associated error bars, is 
shown in Figure 7 as dots. These errors are the statistical ones 
that pertain to this particular dataset. Deeper observations, most 
notably in the Ha F656N filter, would lead to a higher number 
of detections (see the discussion in section 4.2). 

The number of excess stars in each panel of Figure 8 is re- 
ported in Table 4, together with the ratio of the total number 
of stars to that of the excess stars, from where one can see that 
roughly 40% of the 781 Ha emitters have quite good photom- 
etry, i.e. an error smaller than 0.2 mag. This generally implies 
errors in both bands less than 0.14, i.e. a 7 sigma detection. 

It is clear from the discussion above that the objects with 
Ha excess can be divided in two groups, according the the 
strength of the (mF675w - wf656n) excess. If it is so strong that 
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Error range 


Number of 


Ratio 




excess stars 


(total/excess) 


0.05-0.10 


7 + 3 


135 


0.10-0.15 


74 ±9 


16 


0.15-0.20 


205 ± 15 


5 


0.20 - 0.25 


264 ± 18 


4 


0.25 - 0.30 


231 ± 15 


4 



Table 4. Distribution of Ha excess stars in different bins of uncertainty 
on the Ha color index mpevsw - MpeseN (see text and Figure 8). 



no "normal" star is expected to match it, i.e. if the dashed his- 
togram in Figure 7 is 0, then the stars can be pinpointed one 
by one. As it can be seen in Figure 7 the condition that the ex- 
pected histogram be corresponds to (wipevsw ^ 'wpeseN) ?J 0.5 
(orEW(HQf) > 17 A, see appendix A). There are 366 such stars 
out of the total 78 1 with Ha excess (47%). If, on the other hand, 
the expected number of stars (dashed histogram in Figure 7) is 
not zero (i.e. mp675w - 'wpeseN ^ 0-5), excess stars can only be 
identified statistically, but not on a one-to-one basis. 

The distribution of the Ha equivalent widths (EW) for the 
781 stars with statistically significant excess is displayed in 
Figure 9. The observed colour excesses (otf675w -'WpeseN) were 
converted into to equivalent widths using the relation derived in 
appendix A and plotted in Figure A. 1 . The sharp drop at low 
values of the equivalent width is not real, but, rather, an artifact 
of our selection criteria that can only identify emission line ob- 
jects above imp(,75w - wf656n) ^ 0.2, i.e. EW(Ha) > 6 A, the 
weaker ones being lost in the noise. 




20 40 
EW(Ha) [A] 

Fig. 9. Distribution of Ha equivalent widtlis for the 781 stars witli sta- 
tistically significant excess. The observed colour excesses (»Jf675w - 
™f656n) were converted to equivalent widths using the relation de- 
rived in appendix A and displayed in Figure A.l. The dashed verti- 
cal line marks the excess above which we can identify the emission- 
line stars unambiguously (mF675w - »tF656N > 0.5, corresponding to 
EW(Ha) > 17 A). 



3.3.2. Balmer continuum excess 

Objects in the South field exhibit yet another spectral pecu- 
liarity in addition to Ha emission: excess continuum emis- 
sion blueward of the Balmer jump at 3646 A, i.e. between the 
F300W and F450W WFPC2 filters. 

The occurrence of this excess emission can be readily seen 
in Figure 10 where we plot the dereddened (mF3oow,o-'»F450w,o) 
vs (mF45ow,o-'"F8i4w,o) colours for stars with overall good pho- 
tometry (^4 < 0.1, see equation 1). For reference, the expected 
locus from the stellar atmosphere models of Bessell et al (1998) 
forZ = 0.3 -20 and log(g) = 4.5 is shown as a solid line. As we 
have already noticed in section 3.1, there is a general very good 
agreement between the locus of the dereddened stars and the 
expectations from the models atmospheres, with a scatter con- 
sistent with observational errors. Thus, the asymmetric broad- 
ening compared to the models at (mF45ow,o - 'mfsmw.o) - 0.7 
cannot be due to photometric errors (in fact, the broadening to- 
wards redder (mpsoow.o - 'nF450w,o) colours is consistent with 
the errors, the one to the blue is not). Rather, it has to be as- 
cribed to genuine deviations in the objects' spectra relative to a 
bare photosphere, as represented by the solid line. 

The excess Balmer continuum emission measured by the 
(wfsoow.o - 'Wf45ow,o) colour positively correlates with the Ha 
excess, as shown in Figure 1 1 where we plot these two quan- 
tities against each another. The high statistical significance of 
this correlation is confirmed by performing Spearman's test on 
the data (e.g., Conover 1980). The value of the correlation co- 



efficient p = -4.5 implies a probability of less than 7 • 10 that 

the two variables are uncorrelated. The fact that p is negative 
means that (otf675w - »Jf656n) increases as (mF300w,o-«iF450w,o) 
decreases towards more negative values, i.e. the equivalent 
width of the Ha emission increases together with the Balmer 
continuum excess. Whereas the quality of the data does not al- 
low to derive the shape of the correlation, its existence is proven 
with a high statistical significance. 

The fact that the Ha and Balmer continuum excesses cor- 
relate hints at a common origin of the two phenomena. 

3.4. Are the peculiar objects Pre-Main Sequence 
stars? 

To recapitulate, in sections 3.3.1 and 3.3.2 we have described 
the detection in the South field of a population of stars that 
exhibit Ha and Balmer continuum excesses. Also, the fact that 
they correlate with each other points towards a common origin 
for them. The question, then, arises as to what they are and 
what causes their observed spectral peculiarities. 

According to the histogram in Figure 9, these stars have 
Ha in emission with equivalent widths in excess of 6 A. This, 
then, excludes any significant contamination of the sample 
by stars with chromospheric activity, whose equivalent widths 
are smaller than 3 A (e.g. Frasca & Catalano 1994, White & 
Basri 2003). 
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Fig. 10. Dereddened (mpsoow.o - 'Mf450w,o) vs (mF45ow,o - mFsww.o) 
colour-colour diagram for the stars in the South field with a mean error 
in the 4 wide bands ((^4 , sec equation 1 ) smaller than 0. 1 . The theoreti- 
cal locus from the stellar atmosphere models of Bessell et al (1998) is 
shown as a solid line and the expected (wfsoow.o - mF45ow,o) colour 
at (otf45ow,o - WFgMw.o) = 0.7 is marked by the thin vertical line. 
The horizontal dashed lines highlight the region used to select stars 
in Figure 11. 

Instead, we propose that the peculiar objects described in 
section 3.3 are Pre-Main Sequence stars. This is based on the 
following considerations: 

1. As we have seen in section 3.2, there are very young, very 
massive stars in the South field and it is only natural to ex- 
pect the presence of equally young, but lower mass stars. In 
fact, Pre-Main Sequence evolutionary models indicate that 
stars less massive than about 2.5 Mq are still approaching 
the Main Sequence after 5 Myr of evolution and occupy 
the same region of the HR diagram as the field population 
subgiants, which are evolving off the Main Sequence after 
several hundred millions or billions of years. 

2. A class of Galactic Pre-Main Sequence stars, the Classical 
TTauri stars, are observed to have Hot emission with 
equivalent widths up to several hundred Angstroms (e.g. 
Fernandez et al 1995). Moreover, both the Ha and Balmer 
continuum excesses in Pre-Main Sequence stars are thought 
to be linked to the mass accretion process from the circum- 
stellar disk (e.g. Calvet et al 2002). As such, a correlation 
like the one we find between the two quantities, and which 
is illustrated in Figure 1 1, is to be expected. 

We refer the reader to Panagia et al (2000) and Romaniello 
et al (2004) for an analogous population in the field of 
SN 1987A, which is less than 10' away from the South field 
we consider here (see Figure 1). 



Fig. 11. Hot excess vs Balmer continuum for stars in the South field 
with 0.4 < (/Mf45ow,o - '"fsmw.o) < 1 (black dots, see Figure 10). The 
thin vertical line marks the expected (raH3oow.o-'"F45()\v.()) colour in this 
(mF45ow,o - Wfsmw.o) interval from the models of Bessell et al 1998. 
The filled squares represent the median (/«f675w - WFeseN) value in 
the (mF3oow,o - 'Wf45ow,o) bin marked by the rectangles, whose vertical 
extent includes 66% of the stars in each bin. The correlation between 
these two quantities is apparent and is confirmed by the high value of 
Spearman's coefficient p = -4.5 (see text). 



The positions in the HR diagram of the 366 Ha excess stars 
that can be identified unambiguously ((wf675w -wpeseN) > 0.5) 
are shown in Figure 12 together with the evolutionary tracks of 
Siess et al (1997). Luminosity and temperature for these stars 
are computed excluding the magnitude in the F300W filter, 
so as to avoid possible contamination from non-photospheric 
emission'. For comparison, there are roughly 7,500 stars in 
the same area of the HR diagram. As we will discuss in sec- 
tion 3.5.2 some of them might be Pre-Main Sequence stars with 
an Ha emission that is too weak to be detected in the available 
images. Most of them, however, belong to the much older LMC 
population, and the 781 Pre-Main Sequence stars we can iden- 
tify in a statistical sense are outnumbered by a factor of about 
10. 

The comparison with Pre-Main Sequence isochrones by 
Siess et al (1997) shown in Figure 12 confirms the presence of 
a generation of young stars also among the low mass stars (the 
evolutionary tracks indicate masses between 0.8 and 2 Mq for 
them). The brightest and coolest among them could be as young 



' For about 75% of the Ho' emission stars the effective temperatures 
determined with or without the F300W band differ by less than 5%. 
For the remaining 25% of the objects, though, the mean difference in 
Teff is 20%, in the sense that the fit including the F300W band returns 
higher values, with a maximum discrepancy of about 40% for the stars 
with the strongest Balmer continuum excesses. 
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Fig. 12. Hess diagram for the general population in the South field 

(greyscale) and the the 366 stars with unambiguous Ha emission 
in the South field ((mpeTsw - "iveset-i) > 0.5, contours and diamond 
symbols to indicate the position of individual emitters). Luminosity 
and temperature for the emission-line objects are computed exclud- 
ing the magnitude in the F300W filter, as it can be contaminated by 
non-photospheric emission. Pre-Main Sequence isochrones for vari- 
ous ages as computed by Siess et al (1997) and the typical errorbar at 
the luminosity of the emission-line stars are also shown. 



as 1-2 Myr, in agreement with what the upper Main Sequence 
indicates. 

However, the majority of the emission objects seem to have 
ages in excess of 10 Myr. A fraction of the emission objects 
even falls below the locus of the Zero Age Main Sequence, a 
fact that, if taken at face value, would be incompatible with 
their suggested Pre-Main Sequence nature. Nonetheless, let us 
note that: (a) the observed spread in the HR diagram is con- 
sistent with being due to photometric and dereddening errors 
(cfr the errorbar in Figure 12), and (b) the stars with Ha ex- 
cess are, on average, colder than the stars at comparable lu- 
minosity, as expected for Pre-Main Sequence objects. This is 
shown in Figure 13, where the temperature distribution for the 
candidate Pre-Main Sequence stars (solid line) is compared to 
the one of all of the stars (dashed line) in two intervals of lu- 
minosity (see Figure 12): -0.2 < log(L/L0) < 0.15 (panel 
(a)) and 0.15 < log(L/Lo) < 0.5 (panel (b)). Applying the 
Kolmogorov-Smimov (KS) test on the unbiimed data and the 
statistics to the binned data confirm that the distribution are 
conclusively different. In fact, the KS probability that the Ha 
and non-Hor emitters are drawn from the same parent distribu- 
tion is less than 1.8% for the brighter sample (panel (a)) and 
less than 0.04% for the fainter one (panel (b)). Similarly, the 
deviations of the biimed data correspond to values of 16.7 
and 36.5, respectively. 
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Fig. 13. Temperature distribution of the stars for which the Ha emis- 
sion can be identified without ambiguity (solid line) compared to 
the one of all of the stars (dashed line) in two intervals of lumi- 
nosity (see Figure 12): -0.2 < log(L/L0) < 0.15 (panel (a)) and 
0.15 < log(L/Lo) < 0.5 (panel (b)). The Kolmogorov-Smimov test 
on the unbinned data confirms that the distributions in each panel are 
different, with the one for the Ha emitters skewed towards lower tem- 
perature, as expected for Pre-Main Sequence stars. 



Taken at face value, the old ages inferred for the candidate 
Pre-Main Sequence stars are at odds with the age of the bright- 
est stars on the upper Main Sequence stars (~ 5 Myr, see sec- 
tion 3.2). However, at this stage, we cannot claim that this age 
difference is real and that star formation in this field has pro- 
ceeded on different timescale for different mass ranges. In fact, 
a shift of about 10% in temperature, in the sense of making 
the candidate Pre-Main Sequence stars colder, is all that it is 
needed to reconcile the two age determinations. Such a shift 
cannot be ruled out on accounts of the uncertainties in deter- 
mining the stellar parameters from optical bands alone (see the 
discussion in Romaniello et al 2002). Let us also note that a 
younger age for the candidate Pre-Main Sequence stars would 
also help to reconcile them with the current understanding of 
Galactic star-forming regions, where evidences of accretion 
tend to disappear in most Pre-Main Sequence stars after about 
10 Myr (e.g. Calvet et al 2005, but see Romaniello et al 2004 
for evidence in the LMC of vigorous accretion at an age of 
about 14 Myr). 

Converting the measured Balmer continuum excesses de- 
scribed in section 3.3.2 to a mass accretion rate following 
the prescription of Gullbring et al (1998, see also Robberto 
et al 2004), one gets values from several times 10"^ up to 
10"' M0 yr"^ . If sustained for significantly longer than 10 Myr, 
accretion rates of this magnitude imply an accreted mass which 
is a large fraction of the total mass of the star and, in turn, a 
large mass of the circumstellar disk. To this end, a younger age 
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for the candidate pre-Main Sequence stars would also reconcile 
our data with the typical disks observed in the Galaxy (a few 
percent of the stellar mass, e.g. Beckwith 1999 and references 
therein). 

In this respect, however, one should also notice that the 
accretion rates inferred here are affected by an obvious se- 
lection effect, in that only the strongest excesses, those above 
the observational threshold, are detected. In addition, accreting 
T Tauri stars are known to exhibit temporal variation in their 
spectral features (e.g. Herbst et al 1994, 2002). At any given 
time, then, only the highest-accreting stars are detected and in- 
tegrating the measured value over the lifetime of the star may 
lead to a gross overestimate of the total accreted mass. 

3.5. The properties of the young population 

Having made the case for the Pre-Main Sequence nature of the 
Ha and Balmer continuum excess stars, in the next two sections 
we will derive global properties of the young population in the 
South field. 
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3.5.1 . Spatial distribution of stars 

Let us now investigate the spatial distribution of the stars in 
the field. Since there is no obvious clustering of stars, we have 
simply divided the field in four regions, coinciding with the 
four chips of the WFPC2. 

As discussed above, in order to study their spatial distribu- 
tion, Pre-Main Sequence stars need to be identified on an indi- 
vidual basis, not just in a statistical sense. Hence, here we will 
consider only those objects with (mF675w - 'Wf656n) > 0.5. The 
spatial distributions of massive and Pre-Main Sequence stars 
belonging to the same generation are shown in Figure 14. The 
actual numbers are reported in Table 5 together with the one of 
Red Clump stars, which are excellent tracers of the older field 
population with ages in excess of a few hundred milUon years 
(e.g. Faulkner & Cannon 1973). 



Chip Massive stars PMS Red Clump 



PC 


2± 


1.4 


6± 


2.4 


25 


±5 


(PC^WF'' 


9.7 ± 


6.8 


28.8 ± 


11.8 


121 


± 24) 


WF2 


24 ± 


4.9 


100 H 


= 10 


63 


±7.9 


WF3 


29 ± 


5.4 


137 ± 


11.7 


61 


±7.8 


WF4 


9± 


3 


123 ± 


11.1 


51 


±7.1 


Total 


64 = 


t8 


366 ± 


19.1 


200 


± 14.1 


(Total^WF" 


19.9 ± 2.5 


114.4 


±6.0 


62.3 


±4.4) 



" Detections in the PC per unit WF area, i.e. 4.8 times larger 
* Total detections per unit WF area, i.e. 3.2 times smaller than the 
entire field. 

Table 5. Spatial distribution of different types of stars in the four 
chips of the South field. The Pre-Main Sequence stars included here 
are only those that can be identified individually, i.e. those with 
(mF675w - mpeseN) > 0.5. For typographical convenience we report 
the Poissonian uncertainty as the square root of the observed number 
of stars in each bin. 



Fig. 14. South Field: spatial distributions of massive (M > 6 Mq; star 
symbol) and Pre-Main Sequence (M < 2 Mq; open squares) stars 
belonging to the same young generation. As in Table 5, the Pre-Main 
Sequence stars included here are only those that can be identified in- 
dividually, i.e. those with (mpeTsw - wpeseN) > 0.5. It appears that the 
density of hot, massive stars is significantly higher in WF2 and WF3 
than elsewhere. 



Most of the massive stars are confined to the WF2 and 
WF3 chips, i.e. in the direction of NGC 2050 (see Figure 1), 
whereas there is a clear deficiency of them in the PC and WF4 
chips. An inspection to Figure 1 shows that the distribution of 
stars is, indeed, very patchy throughout the entire region and 
that the South field happens to be in one of the densest spots. 
Interestingly, the distribution of hot stars does not seem to cor- 
relate with the reddening distribution. The number of massive 
stars is almost the same in the WF2 and WF3 chips, but the 
extinction is appreciably higher in the former than in the latter 
(see Figure 4 and Table 3), as if the massive stars in the WF2 
chip had not yet dissipated the cloud they were born from. As 
a word of caution, however, let us notice that our method to 
determine the interstellar extinction does not allow one to mea- 
sure the total colurmi density of dust along a line of sight, but, 
rather, the one to the star that is used as a beacon. Therefore 
it is possible that the true dust distribution is rather uniform 
across the two chips, but that stars in the "more reddened" one 
are just more embedded in the dust cloud. The fact that the dis- 
persions in E(B - V) are higher where the mean reddening is 
higher seems to favour this scenario (see Table 3). 

Let us note that the number of low mass Pre-Main Sequence 
stars we identify through their strong Hor excess is only a lower 
limit to their true number. Nevertheless, the spatial distribution 
of the detected stars can be taken as representative of the to- 
tal one. The distribution of young low mass stars belonging to 
the same generation as the massive ones is remarkably constant 
throughout the field, with the exception of the small area cov- 
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ered by the PC chip. As a consequence, the ratio of low-mass 
to massive stars is higher where the density of massive stars is 
lower. This is the same trend as Panagia et al (2000) found in 
the field of SN1987A and, as we will see, it is also the case in 
the West field. Our results, then, provide further evidence for 
the existence of an anti-correlation between high and low mass 
young stars, which indicates that star formation processes for 
different ranges of stellar masses are rather different and/or re- 
quire different initial conditions 

The density of Red Clump stars, which trace the field pop- 
ulation with ages in excess of a few hundred million years (e.g. 
Faulkner & Cannon 1973), is constant over the three WF chips 
and it is only marginally higher in the PC chip. Their density 
in the South field is virtually identical to the value measured in 
the West one and around SN1987A (Panagia et al 2000). This 
further confirms that the older LMC population is uniformly 
distributed, at least on the scale of about 100 pc as probed by 
these three fields. 



3.5.2. The mass function and star formation density 

On the basis of their positions on the HR diagram, we can 
estimate fairly accurate masses for the most luminous stars. 
However, the individual masses of lower luminosity stars of 
the young generation are hard to determine because, taken in- 
dividually, their displacements from the MS are comparable to 
the uncertainties of their effective temperatures and luminosi- 
ties. Therefore, in a discussion of the Initial Mass Function 
(IMF), we have to Umit ourselves to comparing the number 
of stars contained over quite large bins of masses. In particu- 
lar, we have 64 stars with 6 < M/Mq < 40 and the 781 stars 
with masses between 0.8 and 2 Mq identified through their Hor 
excess. Their ratio is consistent with a mass function with a 
power-law slope a - -2.37, which is indistinguishable from 
the "benchmark" Salpeter (1955) value of a = -2.35. 

This result, however, is affected by three sources of incom- 
pleteness that may make the true slope of the stellar mass func- 
tion considerably different from the observed one: (1) First, our 
WFPC2 field is at the outskirts of a concentration of massive 
stars (NGC 2050, see Figure 1) and any mass segregation to- 
ward the center of the cluster could reduce the number of mas- 
sive stars as compared to the one of low mass stars, making the 
observed slope steeper than the actual one. (2) Second, we can 
only identify low mass stars showing strong Ha emission that, 
in Galactic star-forming regions, are only a fraction of the total 
number of low mass Pre-Main Sequence stars. Judging from 
the results of the West field (see section 4), this is a function of 
the depth of the Ha images and can lead to an underestimate of 
the true number of Pre-Main Sequence stars by a factor of two. 
As a consequence, we can only give a lower limit to the true 
mass function slope. (3) Third, Pre-Main Sequence stars are 
intrinsically variable (see, for example Bertout 1989) and their 
true number is presumably greater, possibly as much as another 
factor of 2, than the one determined at any given time. Again, 
this makes the observed mass function shallower than the true 
one. For example, if the true number of Pre-Main Sequence 



stars were 2x2 greater than our direct estimate, the slope of 
the mass function would become a = -3. 

The total mass associated to the most recent episode of star 
formation can also be estimated, but it will be affected by an 
even larger uncertainty because, in addition to a possibly in- 
complete count of stars in the range 0.8-2 Mq, we have to ex- 
trapolate down to an unknown lower mass limit. The lowest 
value of the total mass is obtained adopting a relatively high 
value for the lower mass limit (0.5 Mq, Kroupa 2002) and our 
derived value of the IMF of a = -2.37. Thus we obtain a total 
mass of about 3, 100 Mq. Conversely, using a more canoni- 
cal lower mass limit of 0.1 Mq and the high slope estimated 
above, a = -3, the total mass of the young population would 
be much larger: 47, 000 Mq. Considering that these are quite 
extreme values and that the truth is likely to lie in the middle, 
we can conclude that the total mass is likely to be 12, 000 Mq, 
to within a factor of two (Icr). 

Stars of the young generations were formed over a period 
of time of the order of 20 Myr. Therefore, the star formation 
density in this field turns out to be 0.4 M0yr 'kpc^^. This 
value is intermediate between the typical values found in ac- 
tively star-forming spiral disks (up to 0.1 MQyr 'kpc~^) and 
starburst regions (1 MQyr'^hpc'^ or higher, see, e.g. Table 1 
of Kennicutt 1998). We note that, while the star formation ac- 
tivity in the South field is rather high, stiU it is a couple or- 
ders of magnitude lower than in 30 Doradus also in the LMC 
(100 MQyr-^kpc"^ in the central 10 pc, Kennicutt 1998), the 
most luminous complex in the Local Group. 



4. The West field 

The log of the observations of the West field is reported in 
Table 2. Regrettably, the wide band imaging for this field is 
limited to two filters, and, therefore, the discussion will neces- 
sarily be more limited than in the case of the South field. 

In the West field we have detected 7,978 stars down to 
rnpeisyf - 24, 3,416 of which have a mean error in the two 
available broadband filters smaller than 0.1 mag: 



^2= y^^^izwlfksw <o_i (2) 

This condition, in turn, corresponds to mpfevsw ^ 23. The 
Colour-Magnitude diagram for the F547M (V-band like) and 
F675W (R-band Uke) filters is shown in Figure 15. 

A visual inspection to the diagram immediately shows the 
presence of different generations of stars, covering a wide age 
range. In particular, the brightest Main-Sequence stars (about 
10 Mq) hint to an age of 20 Myr or lower for the younger 
population (see below). As in the case of the South field, the 
presence of an old population is revealed by the Red Clump, 
located at mF547M - 19.5, (mF547M - ^Mpevsw) -0.75, and of a 
well developed Red Giant Branch extending up to mF547M - 17 

and (OTF547M - ?Mf675w) - 1. 
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Fig. 15. mF547M vs (mF547M-nJF675w) ColouT-Magnitude diagram of the 
7,978 stars detected in the West field. The black dots are the 3,416 stars 
with a mean photometric error in the two bands (^2, see equation 2) 
smaller than 0. 1 mag. 

4.1. Reddening and ages 

Since the West field was imaged only in two broad bands, we 
cannot recover reddening, effective temperature and luminosi- 
ties for the individual stars as we did for the South field in sec- 
tion 3.1. However, we can still estimate the average redden- 
ing by fitting prominent features in the CMD such as the Main 
Sequence of the young population and the Red Clump. Using 
a distance modulus of 18.57 ± 0.05 (Romaniello et al 2000), a 
mean value of E(B - V) = 0.25 for the whole field provides a 
good fit to the location of both features (see Figure 16). 

In addition. Figure 16 shows that the Main Sequence is 
rather broad, S(mp547u - wf675w) - 0.2 at mF547M - 20, re- 
flecting the combined effect of differential reddening and age 
spread. The location of the most luminous, hence youngest. 
Main Sequence stars indicates an upper Umit of 20 Myr to their 
age (regrettably, the very bluest star at (mF547M - ^Mpevsw) - 
-0.2 is so heavily saturated that its colour is unrehable). 

The old population is composed of a mix of stars with ages 
from several hundred millions to several billions of years. A 
representative isochrone of 5 Gyr is shown in Figure 16. 

4.2. Pre-Main Sequence stars 

Applying the criteria discussed in section 3.4 we have iden- 
tified 761 ±37 stars with a statistically significant Ha emis- 
sion (see Figure 17). Since in the West field there are only 
two broad bands available, we could not recover the effective 
temperature for the stars. Hence, the expected Ha distribution 
was computed using the relation from the Kurucz (1993) mod- 
els between the orily broad-band colour we have, (mF547M,o - 



Fig. 16. Same as Figure 15, but with the isochrones for 20 Myr 
and 5 Gyrs (Z=Zo/3) overplotted to the data. The isochrones have 
been adapted to the WFPC2 bands from the models of Brocato & 
Castellani (1993) and Cassisi et al (1993) with the model atmo- 
spheres by Bessell et al (1998). A distance modulus of 18.57 was used 
(Romaniello et al 2000), yielding an average reddening of E(B - V) = 
0.25. 



wf675w,o), and (mF675w-WF656N). As it can be seen in Figure 17, 
also in the West field stars with (wf675w - »Jf656n) > 0.5 can be 
identified individually as Ha-emitting objects (264, 35% of the 
total of stars with a statistically significant Ho; excess). 

It is worth noting that, even if the West field represents a 
sparser enviroiraient than the South field (1.6 times less stars 
detected in the former than in the latter), the number of Hor- 
excess objects detected is essentially the same in the two fields 
(A^(South) = 781, A^(West) = 761). This is because flie Ha 
exposures of the West field are much deeper than the ones of 
the South field (a factor of 2.4 in exposure time, see Tables 1 
and 2). In fact, if we increase the Ha photometric errors in the 
West field by a factor of '\/2A ^ 1.5 to make them compara- 
ble to the ones in the South field, our procedure would identify 
only 210 candidate Pre-Main Sequence stars. This result pro- 
vides an estimate of the uncertainties inherent in the study of 
an elusive population such as the one of low-mass Pre-Main 
Sequence stars when they are outnumbered by a much older 
field population. 

The spatial distributions of low mass Pre-Main Sequence 
and massive stars in the West field are compared with each 
other in Figure 18 and Table 6 (the Pre-Main Sequence stars 
are the 264 that can be identified on an individual basis, i.e. 

(»lF675W - »If656n) > 0-5). 

Just as in the case of the South field and the region of 
SN1987A (Panagia et al 2000), the ratio between low and and 
high-mass stars of the same generation exhibits substantial spa- 
tial variations (from 44 in the WF4 chip to roughly half of this 
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Fig. 17. As in Figure 7, but for the stars in the West field. Also here, 
the presence of a prominent population of stars with significant Ha 
excess is evident. The value of the excess above which stars can be 
identified individually is marked with a vertical line. 
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Fig. 18. West field: spatial distributions of massive (M > 6 M©; star 
symbol) and Pre-Main Sequence (M < 2 M©; open squares) stars 
belonging to the same young generation. As in Table 6, the Pre-Main 
Sequence stars included here are only those that can be identified in- 
dividually, i.e. those with (mYeisin - mpeseN) > 0.5. It is clear that 
the candidate Pre-Main Sequence stars are mainly concentrated in the 
WF3 and WF4 chips, i.e. in the direction of SN1987A (see Figure 1). 



value in the WF2 and WF3 chips) mostly due to the variations 
of the number of PMS among the various chips. 



Chip 


Massive stars 


PMS 


Red 


Clump 


PC 


0± 1 


10 ±3.2 


8: 


t2.8 


(PC^WF" 


0+1 


48 ± 15.4 


38.7 


± 13.6) 


WF2 


3+1.7 


66 ±8.1 


49 ±7 


WF3 


4±2 


100 ± 10 


54 


±7.3 


WF4 


2+1.4 


88 ± 9.4 


55 


±7.4 


Total 


9±3 


264 ± 16.2 


166 


± 12.9 


(Total^WF" 


2.8 ± 0.9 


82.5 ±5.1 


51.7 


±4.0) 



" Detections in the PC per unit WF area, i.e. 4.8 times larger. 
* Total detections per unit WF area, i.e. 3.2 times smaller than the 
entire field. 

Table 6. Spatial distribution of different types of stars in the four 
chips of the West field. The Pre-Main Sequence stars included here 
are only those that can be identified individually, i.e. those with 
(mF6i5w - ffJFeseN) > 0.5. 



4.3. The mass function and star formation density 

From Table 5 we see that there are 9 stars with masses between 
6 and 10 Mq and, as discussed in the previous section, a total 
of 761 candidate Pre-Main Sequence stars. These numbers are 
consistent with a mass function slope of or ^ -3.1. 

The same discussion presented in section 3.5.2 can be re- 
peated for the West field. However, here incompleteness is 
much less of a problem, both for massive stars because there 
is no sign of a cluster or group of early type stars near this 
field (see Figure 1), nor for low-mass PMS stars because for 
the West field we have much deeper Ha exposures. 

Thus, adopting a slope a = -3.1, and normalising the 
mass function to the number of massive stars present in the 
field we estimate a total mass of recently formed stars between 
~ 2, 300 Mq and 14, 000 Mq, depending on whether we as- 
sume a lower mass cutoff of 0.5 or 0.1 Mq, respectively. 

Taking the geometric mean of the two extremes as a "most 
probable" value of the total mass, Mtot = 5,700 Mq, and a 
formation time interval comparable to the young star ages, i.e. 
~ 20 Myr, the star formation density in the West field turns out 
to be 0.2 Moyr^'kpc"^. 

5. Summary and conclusion 

In this paper we have presented HST-WFPC2 broad and narrow 
band imaging of two fields with recent star formation activ- 
ity in the Tarantula region of the Large Magellanic Cloud. By 
means of the technique pioneered by Romaniello (1998, see 
also Panagia et al 2000 and Romaniello et al 2004), we have 
identified stars stars of 1-2 Mq with a significant Ha emission: 
781 ± 33 in the South field and 761 ± 37 in the West one. The 
Ha emission positively correlates with Balmer continuum ex- 
cess, indicating a common physical origin for the two phenom- 
ena. The measured Ha equivalent widths (up to several tens of 
A; see Figure 9) are far larger than it can be accounted for by 
normal chromospheric activity (e.g. Frasca & Catalano 1994, 
White & Basri 2003). 

Given the presence of massive young stars in both fields 
and the observed correlation between Ha and Balmer contin- 
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uun excesses , we interpret the emission objects as Pre-Main 
Sequence stars. In this respect, then, they are the equivalent of 
the Classical TTauri stars observed in Galactic star-forming re- 
gions. 

The inherent uncertainty in identifying the Pre-Main 
Sequence stars does not allow a precise determination of the 
mass function of the young population. We have discussed the 
different effects that bias its determination and placed a firm 
lower limit to its slope based solely on Ha-excess stars of 
a = -2.37 and -3.1 in the South and West field, respectively. 
The real mass function could be steeper should the selection 
criterion we have adopted miss a significant number of stars. 

Also the determination of the star formation rate associated 
with the young generation of stars is affected by the selection 
criteria applied to identify the candidate Pre-Main Sequence 
stars, with variations of a factor of a few. Our best estimates 
for the star formation density are of ~ 0.4 M0yr~^kpc~^ in the 
South field and ~ 0.2 Moyr 'kpc^^ in the West field. These 
values are intermediate between what is found for actively star- 
forming spiral disks (less than 0.1 Mq yr~^kpc~^) and starbursts 
(1 M0yr"'kpc"^ or more; Kennicutt 1998). 

The relative spatial distribution of equally young stars with 
different masses, however, is not affected by the bias on the 
selection criteria for Pre-Main Sequence stars and provides 
clues to the mechanisms that lead to star formation. For both 
fields the spatial location of the low mass Pre-Main Sequence 
stars does not follow the one of the massive stars of the same 
young age. This further confirms the results found by Panagia 
et al (2000) for the field of SN1987A, also is in the Tarantula 
region (see Figure 1). The almost anti-correlation of spatial 
distributions of high and low mass stars of a coeval genera- 
tion indicates that star formation processes for different ranges 
of stellar masses are rather different and/or require different 
initial conditions. An important corollary of this result is that 
the very concept of an Initial Mass Function seems not to have 
validity in detail, but may rather be the result of a random pro- 
cess, so that it could make sense to talk about an average IMF 
over a suitably large area, in which all different star formation 
processes are concurrently operating. 

Appendix A: The relation between EW(Hor) and 
m(F675W) - m(F656N) 

In order to compute the relation between the colour excess 
m(F675 W) -ot(F656N) and the equivalent width of the Hor line, 
let us begin by defining Tr(A) as the total response of the sys- 
tem in the broad-band filter R (telescope-i-filter+detector) and 
/('I) = fc+fii^) the total flux from the source (continuum-l-line 
emission). If the contiimum is assumed to be flat across the 
broad band, the detected flux in it is: 

fR= { n(A)f{A)dA = 

JR 

= fc f TR{A)dA + TRiWa) f fi{A)dA 
Jr Jr 



this latter quantity can be replaced with its value at the central 
wavelength of the line, TR{Yla). From the definition of equiva- 
lent width (EW), it foUows that: 



fi(A)dA = EWiRa)-fc 



and, hence: 

fR = fc f TRiA)dA + TRi}la)-fc-EWi}la) 
Jr 



(A.1) 



A similar equation holds for the measured flux in the Ha 
narrow-band filter: 



ha=fc f TnamA + TnaOia)-fc-EWiUa) 

J Ha 



(A.2) 



where TuaiA) is the response curve of the narrow band filter. 

Finally, transforming fluxes to magnitudes leads to the re- 
lation we are looking for, although in an implicit form: 



wr - mua = -2.5 log(^//H„) + A(ZP) 



(A.3) 



where A(ZP) is the difference between the instrumental zero 
points in the broad and narrow band. In our case, we have used 
file WFPC2 filters F675W as R and F656N as Ha and A(ZP) = 
4.312 (Whitmore 1995). 

The magnitude difference (otf675w - wif656n) as a function 
of the Ha equivalent width as computed by combining equa- 
tions (A.l), (A.2) and (A.3) is shown in Figure A.l. 




where we have further assumed the line to be much narrower 
than the filter response curve across it so that, in the last term. 



1 2 3 

log(EW(Ha)) [A] 

Fig. A.l. Colour excess mpevsw - 'Mf656n as a function of the Ha equiv- 
alent width computed from equations (A.l), (A.2) and (A.3). 
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